Uncertainties in the Oi//3/3— decay nuclear matrix elements 
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The nuclear matrix elements M^" of the neutrinoless double beta decay of most 

nuclei with known 2i^/3/?-decay rates are systematically evaluated using the Quasiparticle 
Random Phase Approximation (QRPA) and Renormalized QRPA (RQRPA). The exper- 
imental 2i'/3/3-decay rate is used to adjust the most relevant parameter, the strength of 
the particle-particle interaction. With such procedure the M"" values become essentially 
independent on single-particle basis size, the axial vector quenching factor, etc. Theo- 
retical arguments in favor of the adopted way of determining the interaction parameters 
are presented. It is suggested that most of the spread among the published M°'"s can be 
ascribed to the choices of implicit and explicit parameters, inherent to the QRPA method. 
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1 Introduction 

The observation of Qvj3(3 decay would immediately tell us that neutrinos are 
massive Majorana particles (for reviews see [1, 2, 3, 4]; the issues particularly rele- 
vant for the program of decay search are discussed in [5]). But without accu- 
rate calculations of the nuclear matrix elements quantitative conclusions about the 
absolute neutrino masses and mass hierarchies can barely be reached. Despite years 
of effort there is at present a lack of consensus among nuclear theorists how to cor- 
rectly calculate the nuclear matrix elements, and how to estimate their uncertainty 
(see e.g. [4, 6]). Since an overwhelming majority of published calculations is based 
on the Quasiparticle Random Phase Approximation (QRPA) and its modifications, 
it is worthwhile to try to see what causes the sizable spread of the calculated M^" 
values. Does it reflect some fundamental uncertainty, or is it mostly related to differ- 
ent choices of various adjustable parameters? If the latter is true (and we believe it 
is) can one find and justify an optimal choice that largely removes such unphysical 
dependence? 

In our recent papers [7, 8] we have shown that by adjusting the most important 
parameter, the strength of the particle-particle force so that the known rate of the 
2z//3/3-decay is correctly reproduced, the dependence of the calculated Quf3P nuclear 
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matrix elements M^'^ on the things which are not a priori fixed, in particular, the 
number of included single particle states, the different realistic representations of 
the nucleon G-matrix, the axial vector quenching factor etc., is essentially removed. 
The method has systematically bc!c;n applied to calculate the nuclear matrix c;k;- 
ments M^'^ for most of the nuclei with known experimental 2z//3/3-decay rates and 
arguments in favor of the chosen calculation method have been given. 

In this contribution to the MEDEX'05 we briefly review the ideas and the results 
of [7, 8]. 

2 Details of the calculation of OuiSP decay matrix elements 

Provided that a virtual light Major ana neutrino with the effective mass {mpf}), 

N 

{■^PH) = Pei^e^^'nii , {all > 0) , (1) 

i 

is exchanged between the nucleons the half-Hfe of the Oi//3/3 decay is given by 

J- = G'^Eo,Z)\M'"'\'\{m^f,}\' , (2) 

-1 1/2 

where G'^'^{Eo,Z) is an accurately calculable phase-space factor, and M'^'^ is the 
corresponding nuclear matrix clement. Thus, obviously, any uncertainty in M'^'^ 
makes the value of (m/3/3) equally uncertain. 

The elements of the mixing matrix |J7eiP and the mass-squared differences AmP 
can be determined in oscillation experiments. If the existence of the Qv(3p decay 
is proved and the value of T-^j^ is found, combining the knowledge of \Uei\^ and 
Am^, a relatively narrow range of absolute neutrino mass scale can be determined, 
in most situations independently of the Majorana phases ai [3]. 

The nuclear matrix element M''^^ is defined as 

^'"^ = (rl)' - ^ + ^OT + M^li) (3) 

where (|/)) are the wave functions of the ground states of the initial (final) 
nuclei. The explicit forms of the operators Mp', Mq^ and Mj," are given in Ref. [8, 
9]. We note that for qa = 1.25 nuclear matrix clement M''^'^ coincides with Af""^ 
of our previous work [7]. This parameterization is chosen so that we could later 
modify the value of qa and still use the same phase space factor G°'^{Eo, Z) that 
contains g\ = (1.25)^, tabulated e.g. in Ref. [9]. 

In [7, 8] the quasiparticle random phase approximation (QRPA) and its mod- 
ification, the renormalized QRPA (RQRPA), are used to describe the structure 
of the intermediate nuclear states virtually excited in the double beta decay. We 
stress that in the QRPA and RQRPA one can include essentially unlimited set of 
single-particle (s.p.) states, but only a limited subset of configurations (iterations of 
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the particle-hole, respectively two-quasiparticle configurations), in contrast to the 
nuclear shell model. On the other hand, within the QRPA it is not obvious how 
many single particle states one should include. Hence, various authors choose this 
crucial number basically for reasons of convenience. 

It is well known that the residual interaction is an effective interaction depend- 
ing on the size of the s.p. basis. Hence, when the basis is changed, the interaction 
should be modified as well. For each nucleus in question three singlc-particlc bases 
are chosen in Refs. [7, 8] with the smallest set corresponding to Ihuj particle- hole ex- 
citations, and the largest to about Ahw excitations. The s.p. energies are calculated 
with the Coulomb corrected Woods-Saxon potential. The in Rcf. [8] calculations 
have been performed using G-matrix based only on the Bonn-CD nucleon-nucleon 
potential because as it was shown in Ref. [7] a particular choice of the realistic 
residual two-body interaction potential has almost no impact on the finally calcu- 
lated mean value and variance a of M'"", with the overwhelming contribution to 
a coming from the choice of the single-particle basis size. 

In QRPA and RQRPA there arc throe important global parameters renormaliz- 
ing the bare residual interaction. First, the pairing part of the interaction is multi- 
plied by a factor gpair whose magnitude is adjusted, for both protons and neutrons 
separately, such that the pairing gaps for the initial and final nuclei are correctly 
reproduced. This is a standard procedure and it is well-known that within the BCS 
method the strength of the pairing interaction depends on the size of the s.p. basis. 

Second, the particlc-holc interaction block is renormalized by an overall strength 
parameter gph which is typically adjusted by requiring that the energy of the giant 
GT resonance is correctly reproduced. We find that the calculated energy of the 
giant GT state is almost independent of the size of the s.p. basis and is well repro- 
duced with gph « 1. Accordingly, we use gph — 1 throughout, without adjustment. 

Third, an very important strength parameter gpp renormalizes the particle- 
particle interaction (the importance of the particle-particle interaction for the 
decay was recognized first in [10]). The decay rate for both modes of decay 
is well known to depend sensitively on the value of gpp (in = 1+ channel the 
sensitivity originates as a pronounced effect of variation of the degree of the SU(4)- 
symmetry violation by the particle-particle interaction [10, 11]). This property has 
been used in [7, 8] to fix the value of gpp for each of the s.p. bases so that the known 
half-lives of the 21/(3(3 decay are correctly reproduced. Such an adjustment of gpp, 
when applied to all multipoles , has been shown in [7, 8] to remove much of the 
sensitivity to the number of s.p. states, to the NN potential employed, and even 
to whether RQRPA or just simple QRPA methods arc used. This is in contrast to 
typical conclusion made in the recent past (see, e.g., [12, 4]) that the values of M'°^ 
vary substantially depending on all of these things. 

We believe that the 2;/ decay rate is especially suitable for such an adjustment, 
in particular because it involves the same initial and final states as the Oi/ decay. 
Moreover, the QRPA is a method designed to describe collective states as well as 
to obey various smn rules. Both double-beta decay amplitudes, Oz//3/3 and 2z//3/3. 
receive contributions from many intermediate states and using one of them for 
fixing parameters of QRPA seems preferable. 



Czech. J. Phys. 51 (2001) 



A 



3 



3.0 
2.5 
2.0 

V 

1.0 
0.5 
0.0 



"i 



• RQRPA(g^=1.25) 
■ RQRPA(g^=1.0) 

O QRPA (g^=1.25) 
□ QRPA (g^=1.0) 



C) 



"[] 



c 



f 



76_ 82 96„ 100-. 116„, 128^, 130„ 136 

Ge Se Zr Mo Cd Te Te Xe 



Fig. 1. Average nuclear matrix elements (M' and their variance (including the error 
coming from the experimental uncertainty in M^") for both methods and for all considered 
nuclei. For ""^Xe the error bars encompass the whole interval related to the unknown rate 

of the 2i'/3/3 decay. 



It is well known that the calculated Gamow- Teller strength is larger than the 

experimental one. Formally, this quenching could be conveniently accomplished 
by replacing the true value of the axial current coupling constant gA = 1-25 by a 
quenched value qa — 1-0. To see the dependence on the chosen qa value, both values 
of gA have been used (for all multipoles) . The matrix elements M'°^ calculated for 
the three s.p. bases and a fixed gA are relatively close to each other. 

For each nucleus the corresponding average {M'^'^) matrix elements (averaged 
over the three choices of the s.p. space) is evaluated, as well as its variance a. The 
errors induced in (M' by the experimental uncertainties in M^^^ are added to 
the theoretical ones. 

The calculated Oz//3/3 matrix elements are presented in Fig. 1. There the averaged 
nuclear matrix elements for both methods and both choices of gA are shown along 
with their full uncertainties (theoretical plus experimental). One can see that not 
only is the variance substantially less than the average value, but the results of 
QRPA are quite close to the RQRPA values. Furthermore, the ratio of the matrix 
elements calculated with different gA is closer to unity (in most cases they differ 
only by ^20%) than the ratio of the respective gA squared (1.6 in our case). Again, 
such a partial compensation of the ^(^-dependence has its origin in the adopted way 
of fixing gpp to reproduce M^^^ because one needs smaller gpp for smaller gA- 



4 



A 



Czech. J. Phys. 51 (2001) 



One can qualitatively understand why our choHen proccdiu'e stabilizes the M' ^ 
matrix elements as follows: Is/P'^ matrix elements involve only the 1"*" states in the 
intermediate odd-odd nucleus. The contributions of the 1+ multipole for both modes 
of the /3/3 decay depend very sensitively on due to a proximity to the point of 
the phase-transition in the 1+ channel (corresponding to the collapse of the QRPA 
equations of motion). On the other hand, virtual excitations of many states of 
different multipolarities in the intermediate odd-odd nucleus contribute to the M'^^ 
matrix element, due to the presence of the neutrino propagator. The multipoles, 
other than 1+, correspond to small amplitudes of the collective motion; there is no 
instability for realistic values of ^ 1.0. Hence, they are much less sensitive to 
the value of g^p. By making sure that the contribution of the 1+ multipole is fixed, 
we therefore stabilize the M'^^ value. The fact that RQRPA essentially removes 
the instability becomes then almost irrelevant thanks to the chosen adjustment of 
9m- 

One can see in Fig. 2 of [8] that all multipoles J"^, with the exception of the 
1+ and various very small entries, contribute with the same sign and show the 
essential stability of the partial contributions against variation of the basis size. This 
suggests that uncertainties in one or few of them will have relatively minor effect. 
It is instructive to see separately the effect of the sometimes neglected short-range 
repulsive nuclcon-nucleon repulsion and of the induced weak nuclcon currcuits. One 
can see in Fig. 6 of [8] that the conclusion of the relative role of different multipoles 
is affected by these terms. For example, in ^°°Mo the 1~ multipole is the strongest 
one when all effects are included, while the 2~ becomes dominant when they are 
neglected. 

3 Uncertainties of the Oi//3/3 decay matrix elements 

Ideally, the chosen nuclear structure method should describe all, or at least 
very many, experimental data and do that without adjustments. That is not the 
case of the QRPA or the RQRPA. The interaction used is an effective interaction, 
and various parameters are adjusted. One of our goals in [8] was to show that a 
majority of differences among various QRPA-like calculations can be understood 
and possible convergence of the QRPA results can be discussed. 

A detailed list of 13 main reasons leading to a spread of the published QRPA and 
RQRPA results is given in [8]. Some nuclei, like ^°'^Mo, exhibit more sensitivity to 
these effects than others, such as ^^Ge or *^Se. That is confirmed by our numerical 
studies. The two very important reasons are: 

1) The choice of ^fpp (usually fixed to reproduce the experimental half-lives either 

of /3-decays or of 2z//3/3-dccay) 

2) Whether the two-nucleon short-range correlations (s.r.c.) are taken into account 
(for realistic Qpp values neglecting s.r.c. would lead to a twofold increase in M^"^. 

A moderate spread (which can be as much as tens of %) can originate from 
the difference choices of the mean field, many-body approximations (RQRPA, SC- 
QRPA etc.), the size of the model space, the residual nucleon-nucleon interaction 
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(schematic zero-range or realistic interactions based on the G- matrix), rcnormal- 
ization of the axial-vector coupling constant qa = 1-0 -j- 1.25 and the higher order 
terms (h.o.t.) of the nucleon current (induced pseudoscalar and weak magnetism, 
about 30% reduction in M^"). 

The differences between many calculations are understandable just from the 
way Qpp was fixed, the considered size of the model space, the inclusion of the s.r.c. 
and other minor effects. 

Based on our analysis, we suggest that it is not appropriate to treat all calculated 
Oz^/3/3-decay matrix elements at the same level, as it is commonly done (see e.g., 
[6, 12]), and to estimate their uncertainty based on their spread. Clearly, when 
some authors do not include effects that should be included (e.g. the s.r.c. or the 
h.o.t. in the nucleon current) their results should be either corrected or neglected. 
Some effects are correlated, like the size of the model space and the renormalization 
of the particle-particle interaction. Again, if those correlations are not taken into 
account, erroneous conclusion might be drawn. In our works [7, 8] we have shown 
that our way of fixing the model parameters removes, or at least greatly reduces, 
the dependence of the final result on most of the effects described above. 

Even after all relevant parameters have been carefully fixed, the QRPA is not 
able always to describe well all relevant weak transitions. In particular, it is some- 
times impossible to describe simultaneously the 2j//3/3 decay rate as well as the f}~ 
and (3~^/EC matrix elements connecting the 1+ ground states of the intermediate 
nucleus with the ground states of the final and initial nuclei (^"^Mo is a well known 
example of this problem, see e.g. [19]). That is an obvious drawback of the QRPA 
method; it is never meant to describe in detail properties of non-collective states. 
But that is less relevant for the description of integral quantities that depend on 
sums over many states. 

The calculations of the Oi^/J/J-decay matrix elements by Civitarese and Suhonen 
[12] deserve more comments. The authors performed them within the approach 
suggested in [14] which employs the nucleon current derived from the quark wave 
functions. It is very important to note that in this approach the two 
nucleon s.r.c. are not taken into account. The h.o.t. of nucleon current were 
studied [14] with the conclusion of their minor role (in contrast to our one). 

The extension of the work [14] for the case when gpp is adjusted to reproduce 
the single /3-decay amplitudes was presented in [16]. In [12] the nuclear matrix 
elements are calculated in the same way as in [16], however, the obtained results 
differ significantly (see Table II of [8]) from each other. For some nuclei the difference 
is as large as a factor of two. There is no discussion of this there or in the later 
Refs. [12, 17]. It is noteworthy that the largest matrix element in [12] is found 
for the Ouf3f3-deca.y of ^^^Xe that disagrees with the results of other authors. The 
reduction of the Oi^/3/3-dccay of the ^•'^Xe is explained by the closed neutron shell 
for this nucleus: a sharper Fermi surface leads to a reduction of this transition. 

Altogether, the matrix elements of [16, 12] are noticeably larger than the present 
ones. Most of that difference; can be attributed to the neglect of the s.r.c. and of 
the h.o.t. of the nucleon weak current in these papers. 

In summary of this section we list the main arguments why we believe that 
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the procedure of adjustment used in our works [7, 8] {gpp from the 2j^/3/3-dccay) 
is preferable to the procedure advocated in [17] {gpp from the single /3-decay) by 
Suhonen and in [18] by Suhonen and Civitarese: 

i) The QRPA is never meant to describe in detail properties of non-collective states 

like the single beta decay of the ground state of a nucleus. Thus, it is more 
relevant to fix the QRPA parameters using an integral quantity like the 2i//3/3- 
decay half-life and not the beta decay of a single state. 

ii) We have showed in [8] that the first 1+ state of the intermediate nucleus is not the 

only one contributing to the 2i^/3/3-dccay amplitude. For instance, in A=76 sys- 
tem many excited intermediate 1+ states give comparable contributions (Fig. 8 
of [8]). Thus, to give preference to the lowest state is not well-justified, the entire 
sum is actually what matters. 

iii) The contribution of the 1+ multipole to the M^"^ and the corresponding M^"^ are 
correlated (Fig. 8 of [8]). Making sure that the M^'' agrees with its experimental 
value constrains the 1+ part of the M^"^ as well. The sensitivity of 1+ multipole 
component of the M^" to gpp far exceeds sensitivity of the other multipoles. 
Hence, making sure that the whole 1+ multipole is correct is crucial. 

iv) The approach of fixing gpp to reproduce the /3~ decay of the first 1+ state of 
the intermediate nucleus is limited only for nuclei with A=100, 116 and 128, 
where the groimd state of the corresponding odd-odd nucleus has 1+. At the 
same time, the 2i^/3/3-decay half-life is known now practically for all nuclei of 
experimental interest, i.e., one can use it to fix gpp. 

v) The running sum contributions to the M^" in ''^Ge and ^"^"Mo (Figs. 8,9 of [8]), 
for the different multipoles as well as for the total matrix element do not reveal 
a single state dominance. Thus, it is not obvious that it is best to choose any 
particular state or transition for the adjustment. 

vi) Finally, adjusting gpp to reproduce the 2i^/3/3-decay half-life essentially removes 
the dependence on other parameters. While we have demonstrated that for all 
9 nuclear systems, a similar proof was not given in [17]. 

It clearly follows from the above analysis that the Table 3 of [17] does not reflect 
real physical situation and therefore should be disregarded. 

In the very recent publications of Civitarese and Suhonen [18] serious short- 
comings are claimed in the procedure of fixing gpp adopted by us. Their criticism 
is based on the consideration of negative values of the 2i^/3/3-decay matrix elements 
compatible with the data. But they have overlooked that if gpp is fixed to reproduce 
the single /3-dccay data, the problem of negative values of the corresponding matrix 
elements is present as well. In [8] this criticism of [18] has been refuted by physi- 
cal arguments which explain why solutions corresponding to the negative values of 
2;y/3/3-dccay should not he c;onsidcred. 

It is also worth to mention that in [18] the calculations have been performed in 
the QRPA without consideration of the effect of two-nucleon short-range correla- 
tions. The higher multipolarities are strongly suppressed by the s.r.c. and h.o.t. in 
the nucleon current. Thus there is no more a clear dominance of the 2~ multipolar- 
ity found in [18]. Of course, the inclusion of the s.r.c. and the induced pseudoscalar 
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coupling of nucleoli current in the calculation of the M^'^ is a question of physics 
and not just a matter of taste. 

In addition, a new factor (rUe-R)"^ introduced without any justification by Civ- 
itarese and Suhonen [18] into the definition of the M^'^ will definitely contribute to 
the confusion among the experimentalists working in the field. 

4 Summary and conclusions 

We have shown that the procedure suggested in our recent works, Refs. [7, 8], 
is applicable to essentially all nuclei with known 22^/3/? decay lifetimes. Adjusting 
the strength of the particle-particle neutron-proton force gpp in such a way that 
the experimental 2j^/3/3 decay rate is correctly reproduced removes much of the de- 
pendence on the size of the single-particle basis and whether QRPA or RQRPA 
is used. We have also shown that the quenching of the axial current matrix ele- 
ments, parameterized by the reduction of the coupling constant gA, also leaves the 
resulting 01/(3/3 matrix elements almost unchanged; they become insensitive to the 
variations of parameters describing the short-range nucleon-nucleon correlations as 
well. Thus, the resulting OvPf] matrix elements acquire well defined values, free of 
essentially arbitrary choices. We also present arguments while we believe that the 
chosen procedure of adjusting the interaction is preferable to other proposed ways 
of adjustment. 

The differences in most, albeit not all, published QRPA and RQRPA results 
can be imdcrstood. Comparison between the results of different QRPA/RQRPA 
calculations would be facilitated if authors of future publications specify in detail 
what choices of explicit and implicit adjustable parameters they made, and discuss 
the dependence of their result on their particular choice. By following these sug- 
gestions a consensus among the practitioners of QRPA/RQRPA could be reached 
and most of the spread between the calculated nuclear matrix elements, that causes 
much confusion in the wider physics community, would be shown to be essentially 
irrelevant. To reach a convergence of the results obtained using QRPA/RQRPA 
is clearly just an important step on the way to reliable and correct Ouf3(3 decay 
nuclear matrix elements. 
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